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We study the chiral properties of quenched domain wall fermions with several gauge actions. We demonstrate 
that the nearly translationally invariant modes in the fifth dimension that dominate the residual mass for Wilson 
gauge action can be substantially suppressed using improved gauge actions. In particular, we study the Symanzik 
action, the Iwasaki action, the DBW2 action and compare them to the Wilson action. 



1. INTRODUCTION 



written as 



Domain wall fermions provide a very elegant 
solution to the problem of chiral fermions on the 
lattice; in the limit of infinite length Lg of the fifth 
dimension chiral symmetry becomes exact even 
when the lattice spacing does not vanish. Since 
practical implementations can only deal with fi- 
nite Ls some residual chiral symmetry breaking 
is induced. This breaking goes to zero as Lg goes 
to infinity. It is therefore important to chose the 
parameters of the simulation so that an accept- 
ably small chiral symmetry breaking is achieved 
at a reasonably small Lg. 

Recently, RBC § and CP-PACS |,| reported 
that by using the Iwasaki gauge action, rather 
than the Wilson action, a significant reduction of 
the residual chiral symmetry breaking is achieved 
for the same Lg and lattice spacing a. It is 
instructive therefore to study the effects of the 
gauge action in a systematic way. 

We report here our results from our studies of 
the residual chiral symmetry breaking for the 1- 
loop Symanzik, Iwasaki Q], and DBW2 actions. 
The DBW2 action was introduced in |^ and it 
was shown by QCD-TARO in that it is a good 
approximation of the RG flow on the two dimen- 
sional plane of the plaquette and rectangle cou- 
plings. 

All the above gauge actions can in general be 
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where P^i/ is the standard plaquette in the /i, v 
plane, and i?^^ and C^^a denote the real part 
of the trace of the ordered product of SU(3) link 
matrices along 1x2 rectangles and 1x1x1 paths, 
respectively. For the 1-loop Symanzik action the 
coefficients Cq, Ci, and C2 are computed in tadpole 
improved one loop perturbation theory [Q . While 
C2 = and cq = 1 — 8ci for both Iwasaki and 
DBW2, ci = -0.331 and -1.4067 respectively. 

As a measure of the chiral symmetry breaking 
we use the so called residual mass. The residual 
mass is defined as 
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where is the mid-point chiral symmetry break- 
ing term which appears in the axial Ward identity 
of domain wall fermions, and tmin is sufficiently 
large to avoid short-distance lattice artifacts. 

In addition, for each configuration we have de- 
fined the ratio 
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which is very useful in monitoring the chiral sym- 
metry breaking on each configuration. 
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Table 1 



The simulation parameters for all the actions 



Action 


P 


Mp 


Ms 


Ls 


Wilson II 
Symanzik 
Iwasaki 
DBW2 


6.00 
8.40 
2.60 
1.04 


0.404(8) 
0.411(14) 
0.415(13) 
0.399(11) 


1.8 
1.8 
1.8 
1.7 


12-24 

16 

16 

8-16 



2. MEASUREMENTS - RESULTS 

We measured the light hadron spectrum, the 
residual mass, and the ratio R{t) of Eq. ^ on 
16^^ X 32 lattices. We used 50 Symanzik and 
Iwasaki lattices, and 90 DBW2 lattices. For all 
the actions tested the lattice spacing was tuned 
to ^ 2GeV using the p mass at the chiral 
limit to set the scale. This matches the Wilson 
gauge action at (3 — 6.0. The corresponding (3 
values are given in Table |l|. We have optimized 
the value of the domain wall height M^. For all 
actions except the DBW2, the optimum value of 
Ms was found to be 1.8. For the DBW2 action 
the optimum M5 is 1.7, which is what we used for 
the results reported below. 

In Figure |^ we plot the the ratio of eq. ^ 
summed over time slices 4 to 16. As we can see, 
the Wilson and Symanzik action data (squares 
and stars) show very large fluctuations. But the 
Iwasaki and DBW2 show a baseline of small fluc- 
tuations interrupted by spikes. Presumably the 
large fluctuations observed for the Wilson and 
Symanzik action are due to a very high frequency 
of spikes. The important feature of Figure |l| is 
that the number of spikes is drastically reduced 
for the DBW2 action. It is expected that these 
spikes are spatially localized and occur on config- 
urations that have eigenvalues of the transfer ma- 
trix in the 5th dimension very close to 1. We have 
confirmed these expectations by measuring both 
the spacial structure of the spikes and studying 
the spectral flow of the 4D Wilson Dirac operator 
on these configurations (see talk by T. Izubuchi 
at this conference). 

In Fig. H we show the Lg dependence of the 
residual mass. All data are for a^^ = 2GeV and 
bare quark mass niq = 0.02. The value of is 
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Figure 1. R{t) for bare quark mass 0.02 summed 
over t from 4 to 16. The fancy squares correspond 
to Wilson action, the bursts to Symanzik, the di- 
amonds to Iwasaki and the octagons to DBW2. 



within 10% of .55. Since we are interested in dif- 
ferences by factors of ten, and since we know that 
the nires depends mildly on the quark mass |^J^ , 
the data can be directly compared as is. For com- 
pleteness we also present the CP-PACS residual 
mass results We have confirmed their re- 

sults for = 16 (diamond in Fig. For the 
Symanzik action we ran at Lg = 16 only; the 
residual mass turns out to be smaller by a fac- 
tor of three than the Wilson residual mass at the 
same Lg, and it is larger by a factor of two than 
that of the Iwasaki action. The DBW2 action not 
only gives the smallest residual mass, but also has 
the steepest decrease as a function of Lg. Since 
we have only three points, we cannot claim that 
we know the asymptotic behavior; it is interest- 
ing to note however, that a simple fit gives us 
TTtj-egis) (7" with q ~ 0.6. Shamir's one loop 
perturbative j|] result is q = 0.5. 

For the DBW2 action, in addition to the resid- 
ual mass tests we have also measured the light 
hadron spectrum and pseudoscalar decay con- 
stants at both 2GeV and 1.3GeV. The resid- 
ual mass at 1.3GeV is about a factor of two 
smaller than the Wilson residual mass at 2GeV 
{Lg = 16). This allows simulations at coarser lat- 
tice spacing and at moderate values of Lg with 
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Figure 2. Lg dependence of rriics at a ^ ^ 2GeV. 
The octagons correspond to the DBW2 action, 
the squares (CP-PACS |]) and diamond (RBC) 
to Iwasaki, the burst to Symanzik, and the fancy 
squares to Wilson. 



minimal chiral symmetry breaking. Our spec- 
trum results show that the already very good scal- 
ing of domain wall fermions with Wilson gauge 
action is preserved with the DBW2 action. For 
a detailed presentation of our spectrum data see 
the talk by Y. Aoki at this conference. 

It has been argued ||l^ that the enhancement 
of chiral symmetry breaking is caused by con- 
figurations which also represent changing topol- 
ogy (i.e. small instantons and/or dislocations). 
Our data (Fig. 0) indicate that both the Iwasaki 
and the DBW2 actions suppress these configura- 
tions. In fact, Iwasaki [Q| has argued that his 
choice of the ci coefficient is on the boundary 
separating the space of actions which have sta- 
ble instantons and those which have unstable in- 
stantons. Consequently, when these improved ac- 
tions are used, the question whether the topology 
changes arises. We have measured the topolog- 
ical charge on all the actions we tested, and we 
have concluded that the topology changes with 



a significantly lower rate for the DBW2 action. 
The problem gets more severe as we get closer 
to the continuum limit. The Iwasaki action also 
shows signs of slow topology changing. We be- 
lieve that this slowdown in the topology change 
can be overcome by using over-relaxed heatbath 
algorithm and by simply running for a long time 
before saving configurations. Given the consider- 
able cost of measuring fermionic observables with 
DWF we can afford to do so. 

3. CONCLUSIONS 

Using the DBW2 action, we have essentially 
eliminated the problem of residual chiral symme- 
try breaking for domain wall fermions at 2GeV 
with Lg = 16. The DBW2 residual mass is about 
an order of magnitude smaller than the Iwasaki 
mjcs, and two orders of magnitude smaller than 
the Wilson TOi-os- Assuming a 440MeV string 
tension, the bare DBW2 residual mass is about 
BOKeV. The observed topology change slowdown 
may be handled by using longer running between 
measurements. Finally, our results suggest that 
the DBW2 action should also be useful for over- 
lap fermions, since it probably makes the approx- 
imation to the sign function converge faster by 
eliminating the small eigenvalues of the Hermi- 
tian Wilson Dirac operator. 
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